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Abstract: We examine super symmetric models with mixed modulus-anomaly mediated 
SUSY breaking (MM-AMSB) soft terms which get comparable contributions to SUSY 
breaking from moduli-mediation and anomaly-mediation. The apparent (mirage) unifica- 
tion of soft SUSY breaking terms at Q = /j, m i T not associated with any physical threshold 
is the hallmark of this scenario. The MM-AMSB structure of soft terms arises in mod- 
els of string compactification with fluxes, where the addition of an anti-brane leads to an 
uplifting potential and a de Sitter universe, as first constructed by Kachru et al. The 
phenomenology mainly depends on the relative strength of moduli- and anomaly-mediated 
SUSY breaking contributions, and on the Higgs and matter field modular weights, which 
are determined by the location of these fields in the extra dimensions. We delineate the 
allowed parameter space for a low and high value of tan /3, for a wide range of modular 
weight choices. We calculate the neutralino relic density and display the WMAP-allowed 
regions. We show the reach of the CERN LHC and of the International Linear Collider. 
We discuss aspects of MM-AMSB models for Tevatron, LHC and ILC searches, muon g — 2 
and b —* sj branching fraction. We also calculate direct and indirect dark matter detection 
rates, and show that almost all WMAP-allowed models should be accessible to a ton-scale 
noble gas detector. Finally, we comment on the potential of colliders to measure the mirage 
unification scale and modular weights in the difficult case where fj, m i T ^> Mqut- 
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1. Introduction 



Superstring theory allows for a consistent merging of gravitational physics with quantum 
mechanics, while also containing within it the possibility of describing the gauge interac- 
tions of the Standard Model. This impressive theoretical framework does not, however, 
allow the extraction of predictions of physical phenomena at experimentally accessible en- 
ergy scales. A major obstruction to obtaining a predictive superstring theory is the lack of 
understanding of how the degeneracy associated with the many flat directions in the space 
of scalar fields (the moduli) is lifted to yield the true ground state, the problem being that 
many couplings necessary for the extraction of "observables quantities" are determined by 
the ground state values of these moduli. 

The discovery of a new class of compactifications, where the extra spatial dimensions 
are curled up to small sizes with fluxes of additional fields trapped along these extra di- 
mensions has been exploited by Kachru et al. (KKLT)JlJ to construct a concrete model 
with a stable, calculable ground state with a positive cosmological constant and broken 
supersymmetry. This toy model is based on type-IIB superstrings including compactifi- 
cation of the extra dimensions to a Calabi-Yau orientifold, with fluxes along these extra 
directions. While the background fluxes serve to stabilize the dilaton and the moduli that 
determine the shape of the compact manifold, it is necessary to invoke a non-perturbative 
mechanism such as gaugino condensation on a D7 brane to stabilize the size of the com- 
pact manifold. Finally, a non-super symmetric anti-brane (D3) - included in order to break 
supersymmmetry - also yields a de Sitter universe as required by observations. The KKLT 
construction, which yields an example of a low energy theory that has no unwanted light 
moduli, broken supersymmetry and a positive cosmological constant, may be viewed as a 
starting point for the program of discovering a string ground state that may lead to the 
(supersymmetric) Standard Model at low energies, and which is consistent with various 
constraints from cosmology. 

These considerations have recently motivated several authors to analyze the structure 
of the soft SUSY breaking (SSB) terms in models based on a generalization of the KKLT 
set-up Q. The key observation is that because of the mass hierarchy 

m-moduli > m 3/2 » m S USY (1-1) 

that develops in these models, the SSB terms receive comparable contributions via both 
modulus (gravity) and anomaly mediation of SUSY breaking Q, with their relative size 



parametrized by one new parameter a. The hierarchy (1.1) that leads to mixed modulus- 
anomaly mediated SUSY breaking seemingly allows the moduli to decay early enough not 
to disrupt Big Bang nucleosynthesis. It has, however, been pointed out that these decays 
of the moduli would inevitably produce an unacceptably large number of gravitinos (or 
other sparticles) Q which would subsequently decay to the LSP, unless the gravitino is 
itself heavier than ~ 100 TeV, or the density of particles is reduced by some mechanism 
such as a period of thermal inflation^]. Here, we will assume that such a mechanism is 
operative and that the observed dark matter, in our case the lightest neutralino, is produced 
thermally upon subsequent reheating. Upon integrating out the heavy dilaton field and the 
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shape moduli, we are left with an effective broken supergravity theory of the observable 
sector fields denoted by Q and the size modulus field T. The Kahler potential depends 
on the location of matter and Higgs superfields in the extra dimensions via their modular 
weights rii = (1) for matter fields located on D7 (D3) branes, or rij = 1/2 for chiral 
multiplets on brane intersections, while the gauge kinetic function f a = T la , where a labels 
the gauge group, is determined by the corresponding location of the gauge supermultiplets, 
since the power l a = 1 (0) for gauge fields on D7 (D3) branes p. 1 

Within the MM-AMSB model, the SSB gaugino mass parameters, trilinear SSB pa- 
rameters and sfermion mass parameters, all renormalized just below the unification scale 
(taken to be Q = Mqut), are given (in the convention used in the event generator Isajet 
@)by, 

M a = M s (l a a + b a g 2 a ) , (1.2) 
A ijk = M s {-a ijk a + 7i + Ij + Ik) , (1-3) 
M 2 S ( Ci a 2 + 4a& - 7i ) , (1.4) 



where M s = j^t, b a are the gauge (5 function coefficients for gauge group a and g a are 
the corresponding gauge couplings. The coefficients that appear in dl.2|) - (|1.4|) are given 



by a = 1 - Hi, a ijk = 3-m-rij - n k and & = >^./, -^r ~ Eo'oSo^l/i)- Finally, 
Dij k are the superpotential Yukawa couplings, is the quadratic Casimir for the & th 
gauge group corresponding to the representation to which the sfermion /j belongs, ji is the 
anomalous dimension and 7, = 8ir 2 gffig • Expressions for the last two quantities involving 



the anomalous dimensions can be found in the Appendix of Ref. [11], whose notation we 
adopt here. 2 

The MM-AMSB model is completely specified by the parameter set, 

m 3/2 ,a, tan/5, sign{n), m, l a . (1.5) 

The mass scale for the MSSM SSB parameters is dictated by M s = . The phenomeno- 
logical parameter a, which could be of either sign, determines the relative contributions 
of anomaly mediation and gravity mediation to the soft terms, and as mentioned above 
I of I ~ 0(1) is the hallmark of this scenario. Non-observation of large flavor changing neu- 
tral currents implies common modular weights of particles with the same gauge quantum 
numbers: within this framework, this suggests a common location for these fields in the 
extra dimensions. Grand Unification implies matter particles within the same GUT mul- 
tiplet have common modular weights, and that the l a are universal. We will assume that 
all l a = I and, for simplicity, a common modular weight for all matter particles, but allow 
a different (common) one for the two Higgs doublets of the MSSM. 



1 More specifically, these modular weights for chiral superfields are obtained for examples with toroidal 
compactifications, possibly with singularities Q. This is not generic as Calabi-Yau compactifications allow 
for more general choices: for instance m = | may also be allowed (see Refs. ||). 

2 We note that a defined in Ref. Q differs from the definition in Ref. jnj that we use here by Q Rcf |j| = 
in(M^/m 3 ) a^~' ^ e or ig ma l KKLT construction yields a Rof || — 1, corresponding to a ~ 4.6 for 7713/2 ~ 
1 TeV, in the notation used in this paper. 
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The universality of the l a leads to the phenomenon of mirage unification^, 11 1 of 
gaugino masses. In mirage unification, the splitting of the gaugino masses at Q = Mqut 
is proportional to P a {da)/ga, where (3 a are the beta-functions of their associated gauge 
groups. As the gaugino masses run from Mqut to lower energy scales, the RG running 
exactly cancels the GUT scale mass splitting, leading (for a > 0) to unified gaugino masses 
at some intermediate energy scale Q = fi m - 1T / MguTj the scale of unification of gauge 
couplings. Indeed, the observation of gaugino unification at the mirage unification scale, 

/W = M GUT e- 87r2 /^, (1.6) 

is the smoking gun of the MM-AMSB scenario. If a < 0, /j, m i T > MquTi though one 
would have to continue extrapolation still using MSSM RGEs beyond Mqtjt to discover 
this! We will assume hereafter that I / 0, since I = would be distinguished by a gaugino 
mass pattern as in the AMSB framework. While /x m j r determines la, the (unified) value 
of the the gaugino masses extrapolated to Q = /i m j r is M a (/j, m i r ) = M s x (la), and so 
gives the value of M s (and hence ?Ti 3/ / 2 )- We note here that the soft SUSY breaking scalar 
masses of the first two generations also unify at Q = /i m j r , which allows for experimental 
corroboration of mirage unification of gaugino masses. In fact, for cases with n# + 2n m = 2, 
third generation and Higgs SSB terms also unify at /i m ; r ||. 

Phenomenologically attractive features of the MM-AMSB scenario, also referred to 
as the mirage mediation model, are that it provides natural solutions to 1) the negative 
slepton mass squared problem inherent to AMSB models, 2) the SUSY flavor problem, 
plausibly assuming common location for matter fermions in the extra dimensions, and the 
concomitant universality of their modular weights, and 3) the SUSY CP problem, in that 
(up to O(-^j) corrections) there are no physical CP violating phases in \i, gaugino masses 
and A-terms: moreover, there are mechanisms that also yield a real value for b = Bfi^, 

For these reasons, a number of authors have begun exploring the associated collider 
and dark matter phenomenology of these models. Indeed, in Ref. IJ the phenomena 



of mirage unification was noted, while mass spectra were computed in Ref's flq, 11, 12 



13]. In Ref. |[14f |, it was emphasized that bringing the mirage unification scale down 



close to the weak scale ameliorates fine-tuning problems in supersymmetric theories. In 



Ref. [13], mass spectra were computed in the MM-AMSB scenario for two choices of 
modular weights (ra#, n m ) = (0, 0) and (1, ^). Regions of parameter space where the 
neutralino is the lightest supersymmetric particle, and where its predicted relic abundance 
(assuming it is a thermal relic in standard Big Bang cosmology) is in agreement with 
WMAP measurements [15] were mapped out. Allowing that the observed DM may consist 
of more than one component, 0,% h? ~ 0.13. Collider reaches in MM-AMSB parameter 
space were also presented. In Ref. f\M, aspects of LHC detection for MM-AMSB spectra 



where sparticle masses are roughly degenerate were examined. In Ref. [17], it was pointed 
out that measurement of soft SUSY breaking terms at the LHC and ILC could measure 
the matter field and in some cases the Higgs field modular weights, by taking ratios of 



scalar to gaugino masses at // m i r . In Ref. [18], the neutralino relic abundance and direct 



and indirect dark matter detection rates were presented for four modular weight choices. 
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In this paper, we present allowed parameter space regions for a wide range of Higgs and 
matter field modular weights at low and high tan/3 values. We also compute the neutralino 
relic abundance, and note where it is consistent with WMAP measurements. We present 
the approximate reach of the CERN LHC with 100 fb _1 of integrated luminosity, and the 
reach of a ^/s = 0.5 and 1 TeV international linear e + e~ collider (ILC). We organize the 
parameter space discussion and delineate the WMAP allowed regions in Sec. 2. In Sees. 3-6, 
we present various features of MM-AMSB models from a scan over a wide range of modular 
weight choices, confining ourselves to the regions consistent with WMAP measurements. In 
Sec. 3, we present plots of (g — 2) M and the branching fraction for b — > s'y decays. The value 
of {g — 2)^ is found to rule out most of the MM-AMSB parameter space at large negative 
values of a and \x > 0. In Sec. 4, we comment on features of WMAP allowed MM-AMSB 
models relevant for SUSY searches at the Fermilab Tevatron, CERN LHC and the ILC. In 
Sec. 5, we discuss projections for signals via various direct and indirect searches for dark 
matter. We find that almost all WMAP-allowed parameter space should be accessible to 
ton-size direct dark matter search experiments. Signals from indirect searches (especially 
those from high energy gamma rays) are (very) sensitive to the dark matter halo profile 
in our Galaxy. An optimistic choice for this implies that GLAST would be sensitive to 
essentially the entire parameter space, but a different (yet currently viable) choice suggests 
that many models may be below the level of sensitivity for a wide range of . In Sec. 6, 



we extend our earlier discussion [17|, and consider the possibility of determining mirage 
unification and matter modular weights in the case where a < where /i m i r > Mgut- In 
Sec. 7, we present our conclusions. 

2. WMAP- Allowed Parameter Space of MM-AMSB models 

To facilitate calculations within the MM-AMSB framework, we have incorporated it as 
model line number 9 in the Isajet 7.75 event generator |l(|. For a given set of parame- 
ters Ql.5| ), Isajet runs the measured gauge and Yukawa couplings from the weak scale to 
the GUT scale, where the GUT scale is determined by where gauge couplings g\ = <?2- The 
boundary conditions ([L~2j)- (|l.4j ) are imposed at Q = Mgut, and the values of weak scale 
SSB parameters are obtained by numerically solving the complete set of 26 coupled 2-loop 
renormalization group equations. Sparticle mixing matrices are determined by freezing the 
parameters at a scale Q = yJ m i h ' m i R i while non-mixing soft parameters (that determine 
masses) are frozen out at a scale equal to their value [jlS[l . Next, the RG- improved one-loop 
effective potential is minimized at an optimized scale (which accounts for leading two-loop 
terms), allowing the magnitude of fi to be determined, and complete one-loop radiative 
corrections to sparticle and Higgs masses are obtained. The Yukawa couplings are updated 
due to threshold effects. The SSB parameters evolved back to Q = MguTj of course, 
do not match their input boundary values because of the differences in the upward and 
downward evolution. The procedure is iterated until a stable solution is obtained. The 
neutralino relic density Q^^h 2 , BF(b — ► s'y), (g — 2) M , BF(B S — > and the direct 

dark matter detection cross section a{Z\p) are then obtained using the Isatools package 
|2C]. We interface to DarkSUSY[21| to obtain indirect dark matter detection rates. 
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Throughout our analysis, we take all the l a = 1, and assume that all matter fields have 
a common modular weight n m , but allow the Higgs fields to have a different modular weight 
nn not necessarily equal to n m . We organize our discussion of the models by specifying 
the values of (nn, n m ) where each of these is allowed to take the values, 0, 1/2 and 1. 
The phenomenology of models with intermediate values (1/3 and 2/3) of modular weights 
that may be allowed by more general compactifications is presumably bracketed by the 
nine cases that we study. For any choice of tan/3 and the modular weights, the a — m 3 / 2 
plane provides a convenient panorama for an overview of the phenomenology. We begin 
by delineating the regions of this plane that are allowed by theoretical considerations and 
by the constraint on Q^h 2 from WMAP, starting with the examination of the three n m 
cases with nu = 0. 



2.1 nn = cases 

Our first results are presented in Fig. Here we show the allowed parameter space in 
the a vs. m 3 / 2 plane for nu = and a) n m = 0, tan/3 = 10, b) n m = 0, tan/3 = 30, 
c) n m = 1/2, tan/3 = 10, d) n m = 1/2, tan/3 = 30 e), n m = 1, tan/3 = 10 and f) 
n m = 1, tan/3 = 30. We take \i > and m t = 175 GeV throughout, and require points to 
be compatible with approximate sparticle and Higgs mass constraints from LEP2 searches: 
rriff > 103 GeV, m fl > 95 GeV, + > 120 GeV and m h > 110 GeV f2§. 3 
The white unshaded regions do not lead to an acceptable sparticle mass spectrum due to a 
wrong pattern of electroweak symmetry breaking, signalled by tachyonic sfermion masses, a 
negative value of fj? or one of the Higgs squared masses. The turquoise-shaded region leads 
to a top-squark LSP, while the magenta-shaded region leads to a stau LSP: these regions 
would likely give rise to stable colored or charged relics from the Big Bang, and hence are 
also excluded. The blue-dotted region leads to an acceptable sparticle mass spectrum with 
a neutralino LSP, but in this case the calculated relic density Sl^/i 2 > 0.5, in violation 
of WMAP limits. The green-shaded region has a lower relic density, 0.13 < Qgji 2 < 0.5 
but is also excluded. The red-shaded regions have Q^h 2 < 0.13 in accord with WMAP 
and hence are allowed. Frames a) and b) are repeated, but updated, from Ref. |l3[. As 
noted in Ref. [|13| (and also in Ref. |l^] and the updated Ref. 0]), in frame a) the 
allowed red-shaded region for a > occurs because the t\ is quite light, with ~ m z 1 - 



This is due to the large value of the At parameter occurring in Eq. 1.3. A large A 



feeds into the running of the soft SUSY breaking terms mf and m~ via large values of 
Xf = m~ + m? + m 2 Tj + A 2 in the corresponding RGEs, which then accentuates the 

t L t R n u t, o 

impact of the large top-quark Yukawa coupling in driving these to small values. The 
resultant light top squark t\ enhances neutralino annihilation in the early universe via 
top-squark co-annihilation. There is also a red-shaded region around a ~ —2 where the 
weak scale gaugino masses M\ ~ — so that bino-wino co-annihilation (BWCA) acts to 
reduce the neutralino relic density |24|] to the required level. 



3 While the LEP2 SM Higgs mass constraint > 114.4 GeV translates to a corresponding constraint 
on h as long as tua is large, we require a somewhat lower bound owing to an expected 3-4 GeV uncertainty 
on the theory calculation of mi ; for a discussion of how this is affected if tua ~ 150 — 200 GeV, see Ref. [E3| . 
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In frame b), for tan/3 = 30, once again we see a top-squark-coannihilation region 
adjacent to the turquoise-shaded stop-LSP region. But in this case, at large a and low 
m 3/2 ~ 5 TeV, there is also a region of stau-coannihilation. Note in frame b) that the 
WMAP-allowed BWCA region at a ~ —2 has been eliminated. On these plots we also 
denote the approximate reach of the CERN LHC with 100 fb _1 of integrated luminosity. 
The LHC reach has been evaluated in Ref. [25| and found to extend to m g ~ 3.1 TeV in 



the case where rriq ~ m g (as is the case for MM-AMSB models, see Sec. 4). The 100 
fb _1 LHC reach extends up to m 3/ / 2 ~ 60 TeV, essentially covering all the region with 

rn g ~ rriq ~ 3.1 TeV. The reach of a ^/s = 0.5 and 1 TeV ILC is also indicated. The ILC 
reach is determined mainly by the kinematic limit for or f^f^ pair production 

processes p6| , and is somewhat smaller than the reach of the CERN LHC. 

In frame c), we see how the allowed parameter space changes if instead we take matter 
to live on the D3 — D7 brane intersection, for which the modular weight n m = 1/2. A 
value of n m = 1/2 acts to reduce the magnitude of the Ai parameters in Eq. (|l.3| ), and 



to reduce the moduli contribution to the matter scalar masses in Eq. (1.4). As a result of 
the reduced value of A t at Mqut, there is not so strong of an RG push on the top squark 
soft masses, and so we are left with somewhat heavier top squarks in the mass spectrum, 
and no top squark coannihilation region appears. Instead, new allowed regions appear: 
the vertical band near a ~ 6 corresponds to an A- funnel annihilation region p7[| , where 
2m^ ~ rriA- The upper part of the ^4-funnel is somewhat beyond the 100 fb _1 LHC reach. 
In addition, the upper part of the vertical band near a ~ 4 occurs where the value of 
drops, and the Z\ becomes mixed higgsino-bino dark matter (MHDM). This is due in part 
to the relatively low value of IM3I at Mgtjt? which feeds through the RGEs to yield a low 



11 value (see Ref. [28, 29] for a discussion of the low IM3I dark matter model). The MHDM 
has a large annihilation cross section to W + W~ and ZZ, which acts to reduce the relic 
density, as in the focus point region of the mSUGRA model. The broad allowed region at 
low m.3/2 is where ti ~ M\ ~ M2, and so here we actually find mixed bino-wino-higgsino 



dark matter[3C]. In this region, LHC collider events should be rich in b and r jets, since i\ 
is produced in many cascade decays and t\ — ► bW\, while Z2 also gets produced, followed 
by Z2 — > f\f + f\T. A sample point with a mixed bino-wino-higgsino LSP is shown as Point 
1 in Table 1. This point would likely be excluded due to BF(b — ► S7) constraints (see Sec. 
3) and possibly by direct dark matter search constraints (see Sec. 5). The BWCA region 
at a ~ —2 is also present in this frame. 

In frame d), with the same modular weights as c) but with tan/? = 30, the overall 
picture is similar to c) for < a ~ 5, but the Higgs funnel disappears into the stau LSP 
forbidden region. The entire WMAP allowed region is accessible to the LHC. For a < 0, 
the BWCA region present for the tan/3 = 10 case disappears, leaving no WMAP-allowed 
region. 



4 While the LHC reach was calculated in Ref. [^5| for the case of the mSUGRA model, it was found that 
the reach mainly depends on the squark and gluino masses, and not on their particular decay modes, so 
long as there is ample energy release in the squark and gluino cascade decays. For this reason, we expect 
that the reach results of Ref. [p5| will also apply approximately to the case of MM-AMSB models, the 
exception being where the energy released in the primary decays is suppressed for kinematic reasons Jlifl , 
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Table 1: Masses and parameters in GeV units for four cases of the MM-AMSB model. Also shown 
are predictions for low energy observables, together with the cross section for direct detection of 
dark matter. In all cases, we take m t = 175 GeV. 



In frame e), we show results for modular weights tih = 0, n m = 1 with tan/3 = 10. 
Since n m = 1, in Eq. (1.4) the leading contribution proportional to a 2 for matter scalars 
is absent, and so squark and slepton masses are suppressed, and usually we find tachyonic 
or charged LSPs in parameter space. A narrow region around a ~ 5 — 8 survives where 
~ mfn while the Z\ is a nearly pure higgsino state. This region of parameter space 
has tachyonic slepton masses at the GUT scale. The upper portion of this region is in 
accord with WMAP/LEP2 constraints, but is beyond the reach of both the LHC and ILC, 
but may be accessible to ton-size dark matter detectors. A sample point for a = 7.5 and 
m 3/2 = 60 TeV is shown as Point 2 in Table 1. There are also some faint regions that 
survive for negative a values; the WMAP/LEP2-allowed portion of these has sleptons of 
mass ~ 100 GeV, which leads to so-called bulk annihilation of neutralinos through light 
slepton exchange. A sample bulk annihilation point is listed as Point 3 in Table 1. While 
this point is WMAP/LEP2 allowed, it has a large negative contribution to (g — 2)^ (see 
Sec. 3), and is likely excluded. By flipping the sign of fi, the (g — 2)^ contribution flips sign, 
and gives a large positive contribution: see Point 4 of Table 1. In frame f), no parameter 
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space points lead to a viable SUSY spectrum. This is due to a combination of tachyonic 
GUT scale soft masses plus the downward push of a large r Yukawa coupling, such that 
weak scale tachyonic masses are produced all over parameter space. 

2.2 n H = 1/2 cases 

In Fig. ||, we present an overview of the MM-AMSB parameter space for nu = 1/2 and 
various n m possibilities, for tan/? = 10 and 30, with fj, > 0. Frame a) shows the a vs. m 3 / 2 
plane for n m = and tan/3 = 10. The shading and labelling is as in Fig. |. We see once 
again that a region has appeared where the t\ is the LSP, although this region is smaller 
than in Fig. |]a). The reason again is that the At parameter is quite large (reduced only by 
the nu = 1/2 factor in Eq. (ft. .3D ) at the GUT scale, and again this serves to reduce the top 
squark soft masses via RG evolution. The red-shaded region to the right of the turquoise- 
forbidden region again occurs due to top-squark co-annihilation effects. Along the region 
of high m3/2 and a ~ 4, gg will typically be followed by g — » tt\, with t\ — ► cZi, so that 
gluino pair events at the LHC will be characterized by the presence of a pair of top quarks 
plus two soft charm jets and additional E™ ss . The region at a ~ 2.5 is again a region with 
a Z\ that is a mixed bino-wino-higgsino state. There also appears an ^-annihilation funnel 
at a ~ 5, and a small BWCA region at a ~ —2. The picture for tan/3 = 30 is illustrated 
in frame b), where the tau slepton becomes lighter due to the effect of the large tau-lepton 
Yukawa coupling. In this case, the base of the ^4-funnel gets wider due to the addition 
of stau-neutralino co-annihilation. In addition, a region of bulk annihilation through light 
staus opens up at negative a. 

In frame c), we see results for n m = 1/2 and tan/3 = 10. The stop-LSP forbidden 
region has disappeared owing to the increased matter modular weights. Now there exists 
a broad region around a ~ 3 — 4 with a mixed higgsino-bino Z±, while the low 771,3/2 part 
includes an LSP with a significant wino component. The ^4-annihilation funnel is present 
at a ~ 6 — 7, and a BWCA region occurs at a ~ — 2. If we shift to large tan (3 = 30 in 
frame d), then a large excluded stau-LSP region occurs. To the right of this region occurs 
an area of mixed A-funnel/stau co-annihilation, while the region to the left is characterized 
by a higgsino-LSP with some degree of stau co-annihilation. We note here that some part 
of the ^-annihilation funnel extends out beyond the reach of LHC with 100 fb ; in this 
region, nig ~ niq ~ 3.5 TeV. Any WMAP-allowed regions at negative a have disappeared. 

For 77, m — 1 and tan /3 — 10, portrayed in frame e), almost all allowed regions have 
disappeared. There are a few allowed points to the left of the excluded region at positive a, 
plus a small region of bulk annihilation via light sleptons with mass ~ 100 GeV at a < —5. 
If we increase tan (3 to 30 as in frame f), then only a few points persist at positive a, while 
a small band around a ~ — 6 exists for negative a which is characterized by a bino-like 
LSP and stau co- annihilation. 

2.3 riH = 1 cases 

In Fig. ||, we show the parameter space regions for Higgs modular weight nu = 1, and 
various matter modular weights with tan/3 = 10 and 30. In frame a), with n m = and 
tan/3 = 10, it is noteworthy that the stop-LSP region which is present in Figs. |l]a) an d 
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has disappeared, even though the matter modular weight remains the same. In this case, 
just as in the n m = 1/2, uh = case in Fig. [j]c, the large Higgs modular weight reduces 
the value of A^Mqut) enough so that the top squark soft masses are not driven to such 
low values. The large Higgs modular weight also reduces the GUT scale Higgs mass, which 
also reduces the weak scale A mass[|3l|]. Thus, the A funnel region moves out to somewhat 
higher values of a than found in Fig. ||a). There also remains a large region of higgsino 
and mixed higgsino-wino-bino dark matter around a ~ 3. The BWCA region maintains a 
presence at small, negative a values. The CERN LHC can cover all the parameter space 
shown. Moving to tan/3 = 30 in frame b), we find the base of the A-funnel widened 
somewhat compared to frame a) due to the presence of stau co-annihilation. In addition, a 
region has appeared at low 7713/2 and negative a where neutralino bulk annihilation through 
light staus can occur. 

In frame c), we plot parameter space for n# = 1, n m = 1/2, the same choices adopted 
in the non-zero modular weight analysis of Ref. J13| |. Here, we have an ^4-annihilation 
region at large a, a higgsino region at small positive a, and a BWCA region at small 
negative a. In this case, the LHC covers the entire ^4-funnel, unlike Figs. |l]c) and ||c). In 
frame d) for tan /3 = 30, a large stau-LSP excluded region has appeared. In addition, the 
A-funnel has moved to very large a values, and is very broad, owing to an overlap with 
the stau co-annihilation region. 

Moving to uh = 1 and n m = 1 portrayed in frame e) for tan (3 = 10, we see that all 
the a > region is now excluded due to a stau LSP. In this case, the moduli contribution 
to scalar masses only comes from the mixing term in Eq. 1.4, and so the scalar spectrum 
is similar to AMSB, where sleptons have tachyonic masses. For negative a values, the 
AMSB-moduli interference term in Eq. L4 becomes positive, so some spectra are allowed. 
The only WMAP-allowed points occur at very low m 3 / 2 values, and are typified by bulk 
neutralino annihilation through light sleptons. In the case of large tan/3 = 30 shown in 
frame f), again, no allowed regions appear for a > 0, while some points at large negative 
a are allowed where neutralino annihlation occurs through light slepton exchange in the 
t-channel. This region, we will see, turns out to be largely excluded because it yields a 
large, negative value of {g — 2) M , in contradiction with experimental measurements whose 
discussion we now turn to. 



3. (g - 2) M and BF(b 37) in the MM-AMSB model 

3.1 (g-2)„ 

Current measurements of the muon anomalous magnetic moment show an apparent devia- 
tion from SM predictions. Combining QED, electroweak, hadronic (using e + e~ —* hadrons 
to evaluate hadronic loop contributions) and light-by-light contributions, and comparing 
against measurements from E821 [|32| at BNL, a positive deviation in a M = — 2 2 ^ of 

Aa^ = a e * p - al M = 22(10) x 10~ 10 (3.1) 



-9- 



is reported in the Particle Data Book[33], i.e. a 2.2a effect. Within this framework, 
corrections from Standard Model expectations arise from one loop diagrams with either 
Wi — or Zi — fii t 2 in the loop, possibly accounting for the (rather weak) discrepancy. 

In Fig. | we show a^ USY |4|, the SUSY contribution versus vn^ L arising from the 
various modular weight choices in the MM-AMSB model. Only points with Q^h 2 < 0.13 
are plotted in this figure. The various colors/shapes of the plotting symbols correspond 
to different choices of modular weights, as listed on the figure. We see from Fig. ||a) 
and Fig. ||&) that, as expected, the SUSY contributions from models with large m^ L (and 
concomitantly large and mp, R ) give almost no contributions, since the contributions 
are suppressed by the large sparticle masses. However, depending on the value of tan j3, 
models with ra^ L ranging from about 200 GeV to several hundred GeV can give sizable 
contributions, and even accommodate the theory-experiment deviation. For the larger 
tan /3 case, very light sleptons in fact give too large a SUSY contribution to the muon 
magnetic moment. In addition, several sets of models give large negative contributions to 
a^ USY . These are the models with n m = 1 and any choice of n# but with moderate to 
large values of negative a, wherein matter scalar masses are suppressed, and the points are 
WMAP-allowed by bulk annihilation of neutralinos via light t-channel slepton exchange. 
Even though we take fi > in these models, we know that 



SUSY mjj/iMjtan/? , , 

% M 1 ( } 



and since the weak scale gaugino masses M 12 are usually negative in these cases, a negative 
contribution to a^ USY arises. Thus (even allowing for considerable theoretical uncertainties) 
these models with o^ USY ~ —40 are likely ruled out. 

The situation is more clearly illustrated in Fig. ||, where we plot a^ USY versus a for 
a) tan/3 = 10 and b) tan/3 = 30. We see that models with a < —5 give large negative 
contributions to a^ USY . Thus, when a ~ —5, models with fi < should give Aa^ allowed 
regions. Models where agreement with the measured CDM relic density is obtained via 
BWCA (these do not occur for the tan/3 = 30 cases) yield a smaller value of o^ USY , in 
part because the sleptons are heavier, and in part because of the relative sign between Mi 
and Mi that leads to a negative interference between the chargino and neutralino loop 
diagrams. 

3.2 BF(b -> s 7 ) 

The branching fraction BF{b — * s'y) serves as a strong constraint on SUSY models in 
part because, if msusY — Myy, supersymmetric contributions mediated by Wtfj and bH + 



loops are expected to occur at similar rates[35| to the SM contribution, mediated by a tW 



loop. The measured branching fraction- from a combination of CLEO, Belle and BABAR 
experiments^]- is BF(b — > 57) = (3.55 ± 0.26) x 10 -4 , while the latest SM calculations 
find@ BF{b -> 57) = (3.29 ±0.33) x 10" 4 . Thus, any SUSY contribution to BF(b 57) 
seems highly suppressed. 

The results from the MM-AMSB models for various modular weight choices and with 
n Si h 2 < 0.13 are shown in Fig. ||, where we plot BF(b — * 57) vs. m- g for d) tan/3 = 10 and 
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b) tan P = 30. Many of the MM-AMSB models with m g ~ 1000 GeV predict BF(b -> 57) < 
2 x 10~ 4 , and thus are likely ruled out (unless small flavor-changing contributions to squark 
mass matrices are invoked). Since m g is strongly correlated with rriq (see Fig. ^ below), 
these are models with relatively light squarks (and sleptons). Models with a < that 
are in agreement with the measured CDM relic density due to bulk annihilation via very 
light sleptons (and correspondingly light squarks) typically give too high a prediction for 
BF(b — > S7), and are again likely ruled out. In frame b) with tan p = 30, the contributions 
to BF(b — > sj) can be even more anomalous, and a much larger fraction of models with 
m g ~ 1000 GeV are likely ruled out. 



4. Collider searches in the MM-AMSB Model 
4.1 Fermilab Tevatron pp collider 

It is possible to search for gluino and squark pair production at the Fermilab Tevatron pp 
collider by looking for multi-jet+S™ lss signals. While many Tevatron E™ ms + jets searches 
for models with gaugino mass universality have been pre-empted by LEP2 searches for 
chargino pair production, in models with non-universal gaugino masses, gluino and squark 
searches may still be of interest. Indeed, in the so-called "low M3 dark matter model" 
(LM3DM) p8|] , the Tevatron with 5 fb _1 of integrated luminosity was found to be sensitive 
to m g ~ 200 - 350 GeV in models beyond the reach of LEP2[[29|. The MM-AMSB models 
share with LM3DM models the characteristic that M3 is reduced relative to Mi and M2 at 
the GUT scale, since the AMSB contribution to M3 subtracts from the moduli-mediated 
contribution. 

In Fig. [/], we plot rriu R vs. m g for MM-AMSB models with a) tan P = 10 and b) 
tan P = 30, and all modular weight choices, where we require Vt^h 2 < 0.13 and > 103 
GeV. The diagonal dashed line denotes where m„ B = m g . One lesson from this plot is that 
generic first and second generation squark masses are comparable in mass to m g for all 
modular weight choices; the possibility that rriq S> m g - which occurs in the large rriQ 
region of the mSUGRA model- does not occur here. We also find that in some cases (such 
as nil = I, n m = 0) the value of m g ~ rriq can be as low as ~ 300 GeV, and hence may be 
amenable to Tevatron gluino and squark searches, where the current limit for rriq ~ m g is 



325 GeV within the mSUGRA framework g§ 



It is also possible for Tevatron experiments to search for SUSY via pp — > 
3£ + E™ ss + X events p9| . The clean trilepton search is viable when the "spoiler" decay 



modes Z2 — > Z\h and Z2 — > Z\Z are closed f40||, *-e. when m^ 2 — < Mz and m^. We 
plot in Fig. § the 

m z 2 ~ m Zi mass difference versus rn^ for all modular weight choices 
and WMAP/LEP2-allowed models, for a) tan/3 = 10 and b) tan/3 = 30. The region above 
the horizontal dashed line is where the spoiler decay Z2 — ► Z\Z turns on. We see a large 
fraction of MM-AMSB models of all modular weights are below this line, and thus have 
closed spoiler modes, so that Z2 — > Z\ii is likely to have a significant branching fraction 
(an exception occur if interference effects suppress the leptonic branching decay or if the 
decay Z2 — ► rf\ becomes accessible). Depending on BF(Z,2 — ► Z\U), a portion of the lower 
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range of ~ 100 — 200 GeV and rn-z 2 ~ m z 1 ^ may thus be accessible to Tevatron 
clean trilepton searches. 

Finally, in Fig. |9] we plot the values of nih vs. for WMAP/LEP2-allowed MM- 
AMSB models and all modular weight choices. Here, we see that nth ~ 125 — 130 GeV, so 
that Tevatron Higgs searches may see evidence of Wh, Zh production depending on the 
integrated luminosity achieved. In addition, the t\ mass can range down to ~ 100 GeV 
in MM-AMSB models and may be accessible to Tevatron searches|41]. However, in these 



cases, the — mass gap is usually quite small, in which case the t\ decay products 
will be soft, making hadron collider searches more difficult. 

We note here also that in the MM-AMSB models with < mt and is light, 
some of the conditions for successful electroweak baryogenesis are satisfied. However, in 
Ref. [42], it is noted that the soft term < condition also ought to be satisfied. We 



find that in MM-AMSB models m| > and the light t\ arises in part due to the large 
At term; this apparently contradicts the requirements for successful EW baryogenesis in 
these models. Of course, this would be moot unless there is also a significant CP violating 
phase (~ 0.05) either via the 6-term, or via corrections to other SSB parameters. 

4.2 CERN LHC pp collider 

The CERN LHC pp collider will begin data taking in earnest in 2008 running at = 14 
TeV. The search for supersymmetry by Atlas and CMS will initially focus on gluino and 
squark production followed by cascade decays. It has been shown in Ref. |25|] that in the 
case where rriq ~ m g as in MM-AMSB models, LHC experiments should have a reach to 
nig rsj 3 TeV for an integrated luminosity of 100 fb _1 . Except for possibly cases where 
the sparticle masses are all roughly degenerate (where yUmir ~ 1 TeV) so that the visible 
decay products are significantly softer than the expectation in the mSUGRA model |l6| , 
we would expect a similar reach within the MM-AMSB framework. We can see from Fig. 
[7|that most MM-AMSB models should be detectable at CERN LHC. Exceptions occur m 
the cases of nu = 0, n m = 1 models for tan/3 = 10, and nu = 1/2, n m = 1 models for 
tan/3 = 30. In these instances, the thin allowed region to the left of the excluded region 
for a ~ 5 — 7 gives rise to cases with m g ranging up to 4 — 5 TeV, which is well-past the 
LHC reach for any projected integrated luminosity value. 

Fig. |8] shows also that a large fraction of models contain a mass difference -m^ < 
Mz- Thus, decays of Z2S produced mostly via cascade decays of gluinos and squarks, should 
lead to a measurable dilepton mass edge with m(U) bounded by ~ m z 1 umess the 
branching fraction for the decay is strongly suppressed. This mass edge can serve as a 
starting point for gluino and squark cascade decay reconstruction [43]. 



4.3 Linear e + e collider 

Experiments at a linear e + e~ collider will most easily be able to discover charginos, and 
sleptons (including sneutrinos, if these decay visibly) via their pair production. Moreover, 
unless the mass gaps — or m^—m^ are very small, signals from their production 
should be readily visible over Standard Model backgrounds essentially all the way to the 
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kinematic limit for their production. Specialized analyses allow the signal to be extracted 
for smaller values of the sparticle-LSP mass gap [26]. We show in Fig. |l0| the f\ mass 
versus W\ mass for all MM-AMSB models consistent with WMAP/LEP2. A y/s = 500 
GeV machine should have a reach to or rxif x ~ 250 GeV, while a y/s = 1000 GeV 
machine will have a reach to m^? or ~ 500 GeV, as indicated by the dashed reach 
lines. Thus, while the ILC would certainly see the light Higgs boson h in all MM-AMSB 
cases, it would only be able to access sparticles for a fraction of the models plotted. 



5. Dark matter searches 



In addition to collider searches for sparticles and Higgs bosons, searches for dark matter 
relics from the Big Bang, either via direct detection experiments at underground dark mat- 
ter detectors, or indirectly via searches for neutralino annihilation to high energy neutrinos 
in the core of the sun, or neutralino annihilations to high energy gamma rays or antimatter 
in the galactic halo, may provide an independent signal for new physics. The results of 
these DM searches are sensitive to the composition of the LSP, and combining these with 
information from collider experiments may help us to zero in on the underlying SUSY 
model. 



5.1 Direct dark matter searches 

In the case of particle physics models with imparity conserving supersymmetry and a neu- 
tralino LSP, dark matter detectors look for rare neutralino-nucleus collisions in experiments 
located deep underground, where the neutralino is expected to deposit of order ten keV 
of energy in an elastic scattering event. The neutralino-nucleus scattering cross section 
naturally breaks up into spin-independent (SI) and spin-dependent (SD) pieces, where at 
present the greatest experimental sensitivity is to SI cross sections, since in this case the 
neutralino couples coherently to the entire nucleus, so that the cross section increases as 
A 2 , where A is the mass number of the nucleus. To facilitate comparision between dif- 
ferent detector materials, the spin-independent scattering cross section off a single proton 
is usually used as a figure of merit. We use the IsaRes code (part of the Isajet/Isatools 
package p0|, f0|]) to evaluate this associated direct dark matter detection cross section in 
MM-AMSB models. 



In Fig. 11, we plot the expectation for agj(Zip) as a function of for all nine 
sets of modular weights, retaining only those points consistent with constraints from both 
WMAP and LEP2 data for a) tan (5 = 10, and b) tan (5 = 30. Presently, the most stringent 
limit on this cross section comes from the CDMS collaboration|45|1 , and is shown by the 
top- most contours labelled CDMS. Since the neutralino can transfer the maximum energy 
to the nucleus if its mass is ~ the mass of the nuclear target m-r, the sensitivity is the 
greatest when ~ mx- The maximum energy that can be transferred reduces rapidly 
if mjj <C ttit, and saturates if the neutralino is very heavy, accounting for the general 
shape of this contour (as well as of the projected reach contours discussed below). We see 
that, especially for the larger value of tan (3, the CDMS experiment has already excluded 
neutralinos as heavy as 300-400 GeV for some choices of modular weights. CDMS II is now 
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operating, and is expected to reach cross sections at the 10 -8 pb level which will begin 
probing MM-AMSB parameter space with neutralinos as heavy as 500-800 GeV depending 
on tan j3. CDMS plans an upgrade to 7 supertowers located in the Sudbury underground 
site, where a sensitivity to cross sections at the 10 -9 pb level is anticipated ^(|. To go 
further, ton-scale noble liquid detectors will have to be used. As an example, we plot the 



reach of the proposed Warm Argon Project (WARP) 1400 kg detector [47], which will 
endeavour to probe cross sections as low as 10~ 10 pb. We see that if proposed ton-sized 
noble gas detectors do achieve the projected sensitivity, they should be able to probe the 
bulk of MM-AMSB parameter space unless Z\ is rather heavy. Exceptions where light 
neutralinos (m^ ~ 500 GeV ) would escape detection even at these facilities occur when 
the neutralino is essentially bino-like, and the WMAP relic density is attained via squark 
or slepton co-annihilation, or via BWCA. 

5.2 Indirect dark matter searches 

In addition to direct detection experiments discussed above, there are also indirect searches 
for dark matter that rely on the detection of a flux of high energy particles, or of anti- 



particles, produced via the annihilation of neutralinos in our galactic halo [18]. 



The IceCubeJlJj] and ANTARES[50| experiments are designed to detect high energy 
neutrinos from the annihilation of neutralinos gravitationally trapped in the core of the 
Sun[^] . In Fig. we present rates for detection of Vn — * fj, conversions in the Antarctic ice, 
with the muons being detected in the IceCube experiment. We show results for a) tan (3 = 
10, and b) tan (3 = 30. The reference experimental sensitivity that we show is for a flux of 
40 muons per km 2 per year with > 50 GeV. We see that while many sets of modular 
weights lead to a potentially detectable signal for neutralinos as heavy as 1 TeV {e.g. 
nu = l, n m = 0), there are other models where the signal is always below the observable 
level (e.g. nu = 0, n m = 0). This may be understood if we recognize that the signal at 
IceCube will be largely governed by the higgsino content of the neutralino since this controls 
the neutralino-nucleon scattering cross section (essential for trapping the neutralinos in the 
first place) and, along with the wino content, also affects the neutralino annihilation rate. 
Thus, when we can only have a bino LSP (as, for instance, for nu = n m = case) the 
signal will be small. In contrast, regions with mixed higgsino or mixed wino-higgsino-bino 
DM will typically lead to the largest signals at IceCube. 

Neutralinos in our galactic halo may also annihilate into anti-particles or anti-nuclei, 
thus providing new ways for their indirect detection, though signals from neutrinos not 
pointing from the sun are below their background. Positron detection experiments in- 
clude HEAT||], Pamela[|| and AMS-02 g]; antiprotons may be detected by BESS||, 
Pamela, AMS-02; anti-deuterons by BESSf5§, AMS-02 and GAPS@. For positrons and 
anti-protons we evaluate the averaged differential antiparticle flux in a projected energy 
bin centered at a kinetic energy of 20 GeV, where we expect an optimal statistics and 
signal-to-background ratio at space-borne antiparticle detectors fl58|. We take the exper- 
imental sensitivity to be that of the Pamela experiment after three years of data-taking 
as our benchmark. Finally, the average differential anti-deuteron flux has been computed 
in the 0.1 < Tg < 0.25 GeV range, where Tg stands for the antideuteron kinetic energy 
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per nucleon, and compared to the estimated GAPS sensitivity for an ultra-long duration 
balloon-borne experiment [57]. 

Our projection for the sensitivity of Pamela (p and e + ) and GAPS experiments (D) is 
shown in Fig. [l3|. We show these for tan f3 = 10 (left column) and tan j3 = 30 (right column). 
The expected fluxes depend on the (unknown) details of the neutralino distribution in our 
galactic halo, and are shown for the Burkert profile [59] obtained if it is assumed that the 
central cusp of the DM halo seen in numerical simulations is smoothed out by heating of 
colder particles. This less clumpy halo distribution yields lower fluxes for the anti-particles 
so that our projections for the reach of these indirect search experiments may be regarded 
as conservative. 

We see that typically only a small fraction of the models will be within the reach of 
Pamela, with a somewhat larger fraction within the reach of GAPS. Perhaps, more perti- 
nently, most models with Trig ~ 300—400 GeV are projected to be beyond the reach of these 
experiments. We should mention, however, that this conclusion is sensitive to our assump- 
tion of the Burkert halo profile. If instead we assume a scenario where baryon infall causes 
progressive deepening of the gravitational potential well, a considerably clumpier halo dis- 
tribution is obtained. For the resulting Adiabatically Contracted N03 Halo Model p0|, 
we have checked that we get larger fluxes by factors of ~15-20 (anti-protons), ~5-7.5 
(positrons) and ~10-15 (anti-deuterons) so that the reach is correspondingly increased. 

Indirect detection of neutralinos is also possible via the detection of high energy gamma 
rays [61] from produced by neutralino annihilation in the centre of our Galaxy [62]. These 
will also be searched for by the GLAST collaboration [33|. We have evaluated expectations 
for the integrated continuum 7 ray flux above a EL = 1 GeV threshold in the MM-AMSB 
model. Since these are extremely sensitive to the assumed neutralino halo distribution, 
we show in Fig. 14 results for both the Burkert profile (upper frames) and the Adiabati- 
cally Contracted N03 Halo Model (lower frames), again for tan/3 = 10 (left frames) and 
tan/3 = 30 (right frames). The horizontal line marks the expected sensitivity of GLAST, 
1.0xl0~ 10 cm^s" 1 . We see that while it almost appears that the entire MM-AMSB pa- 
rameter space will be accessible to GLAST for the optimistic projection using the ZV03 
profile, the inference using the Burkert profile is completely different. We need to better 
understand the halo profile before any definitive statement can be made, but at least for 
some models, the reach of GLAST even very conservatively extends out to = 500 GeV. 



6. Testing mirage unification for a < and very large /i m ; r 

The distinguishing characteristic of MM-AMSB models is that the SSB gaugino mass pa- 
rameters, when extrapolated to high scales using one loop RGEs are expected to exhibit 
unification at the scale Q = fi m i T given by Eq. (|1.6[). We emphasize that there is no physical 
threshold at this scale and the phenomenon is therefore referred to as mirage unification. 
For positive values of a, M wca k < /U m i r < Mqut, and it is straightforward to test the 
mirage unification of gaugino mass parameters since the weak scale SSB gaugino masses 
can be extracted from measurements at the proposed electron-positron linear collider. By 
evolving the corresponding weak scale sfermion SSB parameters from the weak scale, the 
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expected concomitant mirage unification of first generation scalar masses- truly a smoking 
gun confirmation of the scenario- may also be revealed |jl7| . An example of the evolution of 
gaugino mass parameters evolved from Mz to Mgut> clearly exhibiting mirage unification, 
is illustrated in Fig.|i~5|a), for the MM-AMSB model with njj = 1, n m = 1/2 with a = 6, 
for which Eq. (1.6) gives fi m i T ~ 4 x 10 10 GeV. As discussed in Ref. fli~7|l , the unified value 



of gaugino masses at // m ; r is given by M s , while the values of the first / second generation 
scalar mass parameters extrapolated to /u m i r are given by m 2 (// m i r ) = (1 — m)Mg, so that 
the ratio of scalar-to-gaugino masses at /i m i r yields 



2 

W 



1-th, (6-1) 



which directly measures that matter field modular weights, and hence their location in the 
extra-dimensional geometry. We should mention though that even with ideal experimental 
measurements, we would not expect mirage unification to be perfect because, 

1. mirage unification is a 1-loop phenomenon, and even if the RGEs are evolved at 1- 
loop, the measured physical masses, and hence the extracted SSB parameters at the 
weak scale include "higher loop effects", and 

2. gauge couplings do not unify exactly because of weak and GUT scale threshold cor- 
rections. 

If a < then /j> m [ T > Mqut an d gaugino and scalar mass parameters would have to be 
evolved beyond Mqut using RGEs corresponding to SU (3) x SU (2) x U(l) MSSM evolution 
to reveal mirage unification. For large negative values of a, this is again straightforward, 
as illustrated in Fig. |l5|fr) for Point 3 in Table 1. We see that the gaugino mass parameters 
exhibit mirage unification beyond Q = .Mqut &t A'mir ~ 

10 20 GeV. Matter scalar mass 
parameters can be similarly evolved to /i m j r in order to extract matter modular weights, 
within errors expected from the extrapolation. 

However, if a is small in magnitude, but negative, we face a new problem because 
/^mir becomes so large that one (or more) of the one-loop-evolved gauge couplings diverges 
for Q < // m i r , making the numerical integration of the RGEs from the weak scale to /t m i r 
difficult. We can circumvent this difficulty because it is possible to analytically integrate 
the one loop RGEs for gaugino masses and gauge couplings to obtain 

where Ki is a constant independent of the energy scale Q, and i labels the gauge group. 
Indeed, if gaugino masses and gauge couplings both unify at M^uis this leads to the 
well-known gaugino mass unification relation. From the experimentally determined weak 
scale values of Mj and gi we can obtain Kf, these can then be used together with l/gf(Q) 
obtained using 

(6.3) 



gf(Q) gf(Qo) 8vr 2 \Q 
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to get the values of l/Mi(Q) at any higher scale. 5 In Fig.|l^c) we show the evolution 
of \/Mi{Q) 2 using the procedure just described for a WMAP/LEP2-allowed case with 
(riH,n m ) = (0,0) with a = —1.635, m 3 / 2 = 25 TeV, tan/? = 10 and /j, > (Point 4 from 
Table 1 of Ref. [jij), which should have a value ^ m ; r ~ 2 x 10 37 GeV, to be compared 
with Q ~ 10 39 GeV in the figure. Notice also that by working with Mi(Q)~ l and l/gf(Q), 
we circumvent the divergence that would have been present at Q ~ 10 27 GeV. The bands 
correspond to a hypothetical uncertainty of ±2% on the extraction of the weak scale values 
of Mi and M 2 , and of ±5% on the extraction of M3. 

Turning to the possibility of checking mirage unification for first /second generation 
scalar masses (for which Yukawa interactions can be neglected) when a < 0, we observe 
that 

m 2 (Q) = Z,- 2 £^M 2 (Q) (6.4) 

is a solution to the 1-loop RGEs for scalar mass squared, where Zi are scale-independent 
constants, and C 2 (fi) is the quadratic Casimir for the & th gauge group that we have en- 
countered below Eq. (|L4|). Since Zi are scale-independent, by taking the difference of 
m?(/i) at the two scales, we obtain, 

™ 2 (^mir) = ™ 2 (^weak) + 2 ^ [M 2 (/i wcak ) - M 2 (/X mir )] , (6.5) 

a a 

where m 2 (/x wea k) is the weak scale SSB parameter for the i th sfermion, which we can obtain 
from its measured mass by removing the D-term. 6 Eq. (|6.4| ), or equivalently, Eq. fl6.5|) can 
then be used to obtain m 2 (^ m ; r ) using the values of gaugino mass parameters at the weak 
and mirage unification scales that we have already obtained, and the corresponding modular 
weights can be extracted. Specifically, we first use the weak scale SSB mass parameters 
to obtain Z^ and then use this value of Zi along with the gaugino SSB parameters at 
Q = /^mir determined above, to obtain m 2 (/i m ; r ). 

We have performed a case study to illustrate this procedure, adopting the point chosen 
in Fig [l5|c), except that we perform the analysis for n m = 0, ^ and 1 to study how well 
the matter modular weights can be extracted in this difficult case. Very optimistically, we 
assume that it will eventually be possible to obtain colored sparticle masses to ±5%, that 
selectron masses will be measured with a precision of ±1%, and finally, that we will be 
able to extract weak scale values of M\ and M 2 to ±2%. From these sfermion masses we 
obtain the corresponding weak scale SSB parameters upon removing the .D-terms using 
cos 2/3 = — 1, and then use Eq. ( |6.4[ ) to get the SSB parameters at the mirage unification 
scale, /i m i r obtained from the intersection of M\ and M 2 (as e.g. in Fig. |i"5|c). The results of 

5 Notice that vanishes (i.e. q?(Q) blows up) when » * , — -^V In ( -3- ) =0, and becomes negative 

gf(Q) v »i r/ gf(Qo) 8?H \Qo J ' ° 

for yet larger values of Q. It is these (unphysical) negative values of gi(Q), obtained from the solution of 
the one- loop RGEs, that must be used in evaluating 1/M;(Q). The ratio of two different gaugino mass 
parameters can flip sign precisely because one of the corresponding gi(Q)s becomes negative. 

6 In practice we may not know the value of tan (3 very well. Fortunately, the change in the D-term is 
negligibly small as long as tan /3 ~ 5, and often even for smaller values of tan /3. In our case study described 
below, we use the D-term with cos 2/3 = — 1. 
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flm = 


n ~ 1 




Q 


""—576 


429 +2 ij 

^—631 


QQQ + 222 

654 


exp 

n Q 


[-7.5, 1.2] 


[-6.9, 1.7] 


[-6.2, 2.4] 




280.1^| 


212.8l? 6 


H0.4™ 


exp 
H er. 


[-0.58, -0.29] 


[0.05, 0.27] 


[0.70, 0.84] 


en 


257.6+Jj 






exp 


[-0.26, -0.13] 


[0.34, 0.43] 


[1.01, 1.03] 



Table 2: Values of sfermion mass parameters at Q — fj, m i T obtained from the intersection of 
Afi(^mir) = ^(/itnir) for the case studies with a — —1.635, m 3 / 2 = 25 TeV, tan/3 = 10, fx > 0, 
Tin = 0, and n m — 0, h, and 1. For each sparticle, we show the value of sign [m? (fi m i T ) ] y/ \ rrif (/i m i r ) | • 
The range corresponds to the projected experimental uncertainties discussed in the text. For each 
sparticle, the second line gives the range of the modular weight consistent with the range of the 
corresponding TO 2 (/u m i r ). The negative sign in the value of m(e^) for n m = 1 really refers to the 
negative sign of m 2 (<=:#). 

our computation are shown in Table |2|. The quoted range on the entries corresponds to the 
variation we obtain by varying the weak scale SSB parameters within their assumed errors. 
For each sparticle, on the second line, we show the range of modular weights compatible 
with the extracted values of the SSB mass squared parameters at the mirage unification 
scale. The following seem to be worthy of note. 

• Even with the optimistic error of only ±5% on the extraction of colored sparticle 
masses, we see that squark SSB parameters are very poorly restricted at the mirage 
scale. The reason is clear from Eq. (|3.5|): squarks and gluinos each have masses 
~ 1000 GeV, so that a very large cancellation between the terms on the right hand 
side is needed to get the left hand side to be ~ (100 — 250 GeV) 2 , the value of the 
matter mass parameters at /u m i r . The effect of higher loop terms, not-exact-unification 
of gauge couplings and various threshold corrections spoil the delicate cancellations, 
leading to a very large error on the extraction of squark modular weights in this case. 

• We see that while slepton masses do appear to unify within ~ 10% at least for 
n m = and n m = 1/2 cases in Table |, the ranges of the modular weights that 
we extract using our procedure are shifted systematically to lower values. We have 
traced this shift to the fact that the value of M a (/x m i r ) ~ 232.5 GeV that we obtain 
from Fig. |l~5|c) is too small by about 10%, causing a systematic shift in the extracted 
modular weights. We had already remarked that // m i r obtained from this figure was 
about 50 times larger than its theoretically expected value. For n m = 1, it appears 
that m 2 (en) and m 2 (ei) unify very poorly. Note, however, that m 2 (// m i r ) is expected 
to be very small in this case, so that we must again have a large cancellation between 
the two terms on the right-hand-side of Eq. fl6,5|) for this to occur. As with the 
squarks, systematic effects then play an important part. It is, however, interesting 
to see that mf{n m \ v ) <C M(/i m j r ), so that the extracted matter modular weights are 
close to unity. 
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We conclude that while our simple procedure provides a good indication of mirage unifi- 
cation of gaugino mass parameters, and a qualitative indication for the mirage unification 
of scalar mass parameters, quantitative extraction of the model parameters have signif- 
icant systematic uncertainties in this difficult case of very large fi m i r - We presume that 
this is because of higher loop and finite correction effects that are present in the weak 
scale SSB parameters. In principle, a determination of all weak scale SSB parameters may 
make it possible to extract "one-loop" weak scale gaugino and first /second generation SSB 
parameters, but until this is possible, quantitative extraction of modular weights (par- 
ticularly with a precision to enable a distinction between toroidal and more complicated 
compactifications) may prove difficult if /i m j r 3> Mqut- 

7. Summary and Conclusions 

Mixed modulus-anomaly mediation of SUSY breaking, also referred to as mirage-mediation, 
provides a new, theoretically well-motivated, and phenomenologically viable framework for 
analyses of the implications of SUSY. The MM-AMSB model is completely specified by 
just three parameters along with the discrete set of modular weights for chiral super-fields 
that are determined by their location in the extra spatial dimensions. The naturally heavy 
moduli and gravitinos can be consistent with cosmology, assuming that there is a second 
late period of inflation that dilutes their density. The observed dark matter may then 
comprise principally of neutralinos that are thermally produced upon reheating after this 
epoch. The framework also provides a plausible solution to the SUSY flavor problem, and 
also ameliorates the SUSY CP problem. 

We have examined the phenomenology of these models under the assumption (moti- 
vated by the absence of large flavor changing neutral currents) that matter superfields all 
reside on either a D7 brane, a D3 brane, or on their intersection: using examples with 
toroidal compactifications, the matter modular weights are fixed to be 0, 1 or |, respec- 
tively. Higgs superfields may lie at a different location, giving us 9 discrete choices for set 
(n,H, n m ) of the matter and Higgs modular weights. More complicated configurations that 
lead to other choices for the modular weights (e.g. 2/3) may be possible, but we expect 
that the phenomenology will then "lie in between" the cases that we examine. 

The phenomenological implications of the model differ from that of other models 
because, (1) the combination of modulus and anomaly-mediated contributions to SUSY 
breaking lead to novel patterns of sparticle masses, and (2) the composition of the lightest 
neutralino differs because the weak scale gaugino masses are not in the ratio ~ 1 : 2 : 6 
predicted by the gaugino mass unification condition, as is the case in many models. These 
features of the spectrum make it possible to obtain a neutralino thermal relic density con- 
sistent with WMAP measurements in a variety of interesting ways that have been proposed 
in previous bottom-up studies, not all of which are possible in the much studied mSUGRA 
model. These include co-annihilation with staus or stops, resonance annihilation via A and 
H, bino-wino co-annihilation, mixed higgsino and mixed bino-wino-higgsino DM and low 
IM3I (mixed-higgsino) DM, but not all mechanisms are realized for every choice of matter 
and Higgs modular weights. Interestingly, mixed wino DM (without also a significant hig- 
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gsino content), or of mixed higgsino DM via non-universal Higgs mass parameters, does 
not seem to be realized within this framework. 

Within this framework, it may be possible for Tevatron experiments to discover super- 
symmetry in portions of the parameter space not accessible to LEP2 searches. Possibilities 
include searches via the E™ 1SS channel for say (n#, n m ) = (1, 0), where gluinos and squarks 
as low as about 300 GeV may have escaped detection up to now, searches for light stops 
and, depending the branching ratio for Z2 — > IIZ\ decay, trilepton events from W1Z2 pro- 
duction. As can be seen from Fig. Q-Fig. ||, assuming 100 fb _1 of integrated luminosity, 
LHC experiments should be able to detect SUSY signals over most of the parameter space 
of the model consistent the relic density measurement, with some exceptions for uh = 
or 1/2, and n m = 1. Experiments at a 1 TeV linear collider will be able to make precision 
measurements over part of the region probed by the LHC. Interestingly, these experiments 
will be able to explore charginos and neutralinos if the measured relic density arises due 
to BWCA: this is important because the small gap, — nig , may make exploration of 
charginos difficult at the LHC. 

We have also examined prospects for direct and indirect detection of DM for the nine 
cases of modular weights. Our results for direct detection are shown in Fig. 11, where we 
see that the CDMS experiment has already begun to exclude a portion of the parameter 
space for as heavy as 300 GeV, for the larger value of tan/3. It is also striking that 
proposed ton size noble gas detectors (we use the proposed Warm Argon Project with 
1,400 kg of argon as a benchmark) will be able to explore most of the parameter space 
of this model, even if is as heavy as ~ 1000 GeV. Prospects for indirect detection 
are somewhat less certain. Models with significant higgsino components yield the largest 
signals, while models with a bino LSP will generally not give an observable signal. There 
may be detectable signals from high energy muons from the sun at IceCube, as well as from 
anti-particles at Pamela and Gaps, or from gamma rays in GLAST. There is, however, 
considerable theoretical uncertainty (factor 10-20) in the prediction of the anti-particle 
signals from the unknown clumping of neutralinos in our galactic halo, and up to four orders 
of magnitude uncertainty for signals at GLAST, so that projections for these experiments 
should be viewed with care. Since this clumping provides a common uncertainty for these 
experiments, observing a signal in one of these will greatly facilitate a much more certain 
prediction for the other experiments, assuming that the observed signal is from neutralinos 
annihilating in our galactic halo. 

Finally, we have explored the possibility of testing mirage unification of SSB param- 
eters, which is the smoking gun for this framework. The case for small or intermediate 
mirage unification scale had already been explored and it had been shown that mirage uni- 
fication of both gaugino and scalar mass parameters would be readily testable, and allow 
extraction of the modular weights. Here, we have explored the much more difficult case 
that occurs if a is small in magnitude, but negative, so that // m j r 3> Mgtjt- In this case, we 
find that mirage unification of gaugino masses should still be testable, although extraction 
of underlying values of fjrmw and M s has significant systematic uncertainty due to the large 
"distance" between weak and mirage scales. Combining measurements from the LHC with 
measurements that will be possible with initial runs of a linear collider may allow a quali- 
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tative test of mirage unification of scalar masses, but unlike the case of intermediate values 
of /x m ir, extraction of matter modular weights is fraught with systematic uncertainty. 

In summary, we have made an exhaustive exploration of the phenomenology of the 
MM-AMSB model, a promising new framework for SUSY phenomenology. The novel 
patterns of sparticle masses that are possible allow various mechanisms that lead to a DM 
relic density in agreement with observations to be incorporated into a top-down framework, 
something not previously possible. This can qualitatively change the correlations between 
various experimental signals from expectations in the much studied mSUGRA framework. 
In view of its attractive and economic theoretical structure as well as its attractive phe- 
nomenological attributes, we believe it is at least as attractive (and theoretically, perhaps, 
more compelling) as mSUGRA as a framework for analyses of SUSY. 
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Gravitino mass vs. a, tanp-10, n H -0, n m -1 
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Figure 1: Allowed regions of parameter space of the MM-AMSB model with modular weight 
riH = and a), b) n m = 0, c), d) n m — 1/2 and e), f) n m = 1. Frames a), c) and e) have tan/3 = 10 
while b), d) and f) have tan/3 = 30. The plots also show the expected thermal neutralino relic 
density from the Big Bang, and the reach of the CERN LHC for 100 fb _1 and for a yfs = 500 and 
1000 GeV linear e+e" collider. We take \i > and m t = 175 GeV. 
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Figure 2: Allowed regions of parameter space of the MM-AMSB model with modular weight 
Tin = 1/2 and a), b) n m — 0, c), d) n m — 1/2 and e), f) n m = 1. Frames a), c) and e) have 
tan j3 = 10 while b), d) and f) have tan (3 = 30. The plots also show the expected thermal neutralino 
relic density from the Big Bang, and the reach of the CERN LHC for 100 fb _1 and for a ^fs = 500 
and 1000 GeV linear e + e~ collider. We take fi > and m t = 175 GeV. 
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Figure 3: Allowed regions of parameter space of the MM-AMSB model with modular weight 
Tin = 1 and a), b) n m — 0, c), d) n m — 1/2 and e), f) n m = 1. Frames a), c) and e) have tan/3 = 10 
while b), d) and f) have tan/3 = 30. The plots also indicate the neutralino relic density, and the 
reach of the CERN LHC for 100 ffcr 1 and a y/s = 500 and 1000 GeV linear e+e" collider. We also 
take fi > and mt = 175 GeV. 
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Figure 4: Plot of vs. mp, L in the MM-AMSB model for various modular weight choices 

with [i > and mt = 175 GeV. Only models with Qgji 2 < 0.13 are retained. In frame a), we show 
results for tan/3 = 10, while in frame 6) we take tan/3 = 30. 
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Figure 6: Plot of BF(b — ► sj) vs. nig in the MM-AMSB model for the different modular weight 
choices, for a) tan/3 = 10, and b) tan/3 = 30. Only models with Vt^ h? < 0.13 are included. We 
take /i > and m t = 175 GeV. 
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Figure 7: Plot of ma R vs. m~ g in the MM-AMSB model for various modular weight choices. Only 
models with £1% h 2 < 0.13 and LEP2 constraints satisfied are included. In frame a), we plot models 
for tan/3 = 10, while in frame b) we plot models for tan/3 = 30. 
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Figure 8: Plot of m^ 2 — vs. in the MM-AMSB model for various modular weight 

choices. Only models with O^/i 2 < 0.13 and LEP2 constraints satisfied are included. In frame a), 
we plot models for tan (3 — 10, while in frame &) we plot models for tan (3 = 30. 
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Figure 9: Plot of vs. mi in the MM-AMSB model for various modular weight choices. Only 
models with Vl? h 2 < 0.13 and consistent with LEP2 constraints satisfied are included. In frame 
a), we plot models for tan/3 = 10, while in frame b) we plot models for tan/3 = 30. 
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Figure 10: Plot of vs. in the MM-AMSB model for various modular weight choices. 

Only models with £1% h 2 < 0-13 and consistent with LEP2 constraints satisfied are plotted. In 
frame a), we plot models for tan/3 = 10, while in frame b) we plot models for tan/3 = 30. 




Figure 11: Plot of asi{Z\p) vs. m^^ in the MM-AMSB model for various modular weight choices 
for a) tan (3 — 10, and b) tan /3 = 30. Only models with £1% h 2 < 0.13 and LEP2 constraints satisfied 
are included. The region above the contour labelled CDMS is excluded by the non-observation 
of a signal in the CDMS experiment. The other contours show the projected reach of CDMS 
(labelled CDMS II) along with the projections for its proposed upgrade (superCDMS) and for the 
Warm Argon Project (WARP 1400), taken here as a representative of proposed ton-sized noble gas 
detectors. 
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Figure 12: The integrated muon flux and the projection for the ultimate reach of the IceCubc 
experiment for the various choices of Higgs and matter modular weights as labelled on the figure 
for a) tan/3 = 10, and b) tan/3 = 30 with fj, > and rrit = 175. The points above the horizontal 
dashed line are projected to be within the reach of IccCube, assuming the sensitivity is given by 40 
events/km 2 /yr for > 50 GeV. Only points consistent with constraints from WMAP and LEP2 
are included in the figure. 
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Figure 13: Expected fluxes and projected sensitivities of the Pamela experiment for the detection 
of anti-protons or positrons, and of the GAPS experiment for the detection of anti-deuterons, 
from neutralino annihilation in our galactic halo, assuming the Burkert profile for the halo DM 
distribution for the various choices of Higgs and matter modular weights as labelled on the figure. 
We take /i > and mt = 175 and show results for tan/3 = 10 (left column), and tan/3 = 30 (right 
column). These experiments should be sensitive to points above the horizontal dashed lines with 
criteria described in the text. Only points consistent with constraints from WMAP and LEP2 are 
included in the figure. 
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Figure 14: The integrated flux and projected sensitivity of the GLAST experiment for the 
detection of gamma rays with E 1 > 1 GeV from neutralino annihilation in the centre of our Galaxy, 
for the the Burkert profile (upper frames) and for the Adiabatically Contracted N03 halo profile 
for the galactic DM distribution, for the various choices of Higgs and matter modular weights as 
labelled on the figure. We take fj, > and m t = 175 and show results for tan/9 = 10 (left column), 
and tan/? = 30 (right column). GLAST should be sensitive to points above the horizontal dashed 
lines with criteria described in the text. Only points consistent with constraints from WMAP and 
LEP2 are included in the figure. 
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Figure 15: Plot of evolution of gaugino masses versus energy scale Q from Q — Mz to Q = Mqut 
and beyond for three cases of MM-AMSB model: a) (nH 7 n m ) = (1, ^) with a = 6 and m 3 / 2 = 12 
TeV, 6) {n H ,n m ) = (0, 1) with a = -10, m 3/2 = 4 TeV and c) (n H ,n m ) = (0, 0) with a = -1.635 
and TO3/2 = 25 TeV. All figures take tan/3 = 10 and fi > 0. We adopt error projections of ±2% on 
Mi and M2, and ±5% on M3 for frame c).Note that we continue MSSM evolution even beyond the 
scale Q = A/gut for reasons that we discuss in the text. 
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